Entangled mixed-state generation by twin-photon scattering 
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We report experimental results on mixed-state generation by multiple scattering of polarization- 
entangled photon pairs created from parametric down-conversion. By using a large variety of scat- 
tering optical systems we have experimentally obtained entangled mixed states that lie upon and 
below the Werner curve in the linear entropy-tangle plane. We have also introduced a simple phe- 
nomenological model built on the analogy between classical polarization optics and quantum maps. 
Theoretical predictions from such model are in full agreement with our experimental findings. 
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I. INTRODUCTION 

The study of spatial, temporal and polarization corre- 
lations of light scattered by inhomogcneous and turbid 
media has a history of more than a century Due 
to the high complexity of scattering media only single- 
scattering properties are known at a microscopic level 
. Conversely, for multiple-scattering processes the em- 
phasis is mainly on macroscopic theoretical descriptions 
of the correlation phenomena [|| . In most examples of 
the latter 0, [E IE 0] *he intensity correlations of the in- 
terference pattern generated by multiple-scattered light 
are explained in terms of classical wave- coherence. On 
the other hand, the recent availability of reliable single- 
photon sources has triggered the interest in quantum cor- 
relations of multiple-scattered light 0] . Generally speak- 
ing, quantum correlations of scattered photons depend 
on the quantum state of the light illuminating the sam- 
ple. In Ref. [U , spatial quantum correlations of scattered 
light were analyzed for Fock, coherent and thermal input 
states. 

In this paper we present the first experimental re- 
sults on quantum polarization correlations of scattered 
photon pairs. Specifically, we study the entangle- 
ment content in relation to the polarization purity of 
multiple-scattered twin-photons, initially generated in a 
polarization-entangled state by spontaneous parametric 
down-conversion (SPDC). The initial entanglement of the 
input photon pairs will in general be degraded by mul- 
tiple scattering. This can be understood by noting that 
the scattering process distributes the initial correlations 
of the twin-photons over the many spatial modes excited 
along the propagation in the medium. In the case of 
spatially inhomogeneous media the polarization degrees 
of freedom are coupled to the spatial degrees of freedom 
generating polarization dependent speckle patterns. If 
the spatial correlations of such patterns are averaged out 
by multi-mode detection, the polarization state of the 
scattered photon(s) is reduced to a mixture, and the re- 
sulting polarization-entanglement of the photon pairs is 
degraded with respect to the initial one. A related theo- 
retical background was elaborated in [E [l(| • 

This paper is structured as follows: In section II we 



report our experiments on light scattering with entan- 
gled photons. First, we present our experimental set-up 
and briefly describe the many different optical systems 
that we used as scatterers. Next, we show our experi- 
mental results. The notions of generalized Werner and 
sub- Werner states are introduced to illustrate these re- 
sults. In section III we introduce a simple phenomeno- 
logical model for photon scattering that fully reproduces 
our experimental findings. Finally, in section IV we draw 
our conclusions. 



II. EXPERIMENTS ON LIGHT SCATTERING 
WITH ENTANGLED PHOTONS 

A. Experimental set-up 

Our experimental set-up is shown in Fig.[TJ A Krypton- 
ion laser at 413.1 nm pumps a 1 mm thick (3 — 
BaB204 (BBO) crystal, where polarization-entangled 
photon pairs at wavelength 826.2 nm are created by 
SPDC in degenerate type II phase-matching configura- 
tion Single- mode fibers (SMF) are used as spatial 
filters to assure that each photon of the initial SPDC 
pair travels in a single transverse mode. Spurious bire- 
fringence along the fibers is compensated by suitably ori- 
ented polarization controllers (PC). The total retardation 
introduced by the fibers and walk-off effects at the BBO 
crystal are compensated by compensating crystals (CC: 
0.5 mm thick BBO crystals) and half-wave plates (A/2), 
in both signal and idler paths. In this way the initial two- 
photon state is prepared in the polarization singlet state 
|^ s ) = (\HV) - \VH))/V2, where H and V are labels 
for horizontal and vertical polarizations of the two pho- 
tons, respectively. The experimentally prepared initial 
singlet state p° xp has a fidelity [l2[ with the theoreti- 
cal singlet state p s = \ip s )(ip s \ of F(p s ,pfP) ~ 98%. In 
the second part of the experimental set-up the idler pho- 
ton passes though the scattering device Ta before being 
collimated by a photographic objective (PO) with focal 
distance / = 5 cm. The third and last part of the experi- 
mental set-up, consists of two tomographic analyzers (one 
per photon), each made of a quarter- wave plate (A/4) 
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FIG. 1: (Color online) Experimental scheme: After singlet 
preparation, the idler photon propagates through the scatter- 
ing system Ta ■ The polarization state of the scattered photon- 
pairs is then reconstructed via a quantum tomographic pro- 
cedure (see text for details). 



II scattering systems are made of birefringent media, 
which introduce an optical axis that breaks polarization- 
isotropy. Birefringence can be classified as "material bire- 
fringence" when it is an intrinsic property of the bulk 
medium (for example a birefringent wave-plate), and as 
"topological birefringence" when it is induced by a spe- 
cial geometry of the system that generates polarization- 
anisotropy, an example of a system with topological bire- 
fringence is an array of cylindrical particles. Finally, type 

III scattering systems are made of dichroic media that 
produce polarization-dependent photon absorbtion. Ex- 
amples of such devices are commonly used polarizers. A 
systematic characterization of all the scattering devices 
that we used was given in Ref. 



C. Experimental results 
in the tangle versus linear entropy plane 



followed by a linear polarizer (P). Such analyzers per- 
mit a full tomographic reconstruction, via a maximum- 
likelihood technique [II], of the two-photon state. Ad- 
ditionally, interference filters (IF) put in front of each 
detector (AA = 5 nm) provide for bandwidth selection. 
Detectors Da and Db are "bucket" detectors, that is they 
do not distinguish which spatial mode a photon comes 
from, thus each photon is detected in a mode-insensitive 
way. 



B. Scattering devices 

All the scattering optical systems that we used were 
located in the path of only one of the photons of the 
entangled-pair (the idler one), as shown in Fig. [TJ For 
this reason, we refer to such systems as local scatterers. 
Such scatterers can be grouped in three general categories 
accordin g to the optical properties of the media they are 
made of jlH]: 

Type I Purely depolarizing media, or diffusers. Such 
media do not affect directly the polarization state 
of the impinging light but change the spatial dis- 
tribution of the impinging electromagnetic field. 

Type II Birefringent media, or retarders. These media 
introduce a polarization-dependent delay between 
different components of the electromagnetic field. 

Type III Dichroic media, or diattenuators. Such me- 
dia introduce polarization-dependent losses for the 
different components of the electromagnetic field. 

Type I scattering systems produce an isotropic spread 
in the momentum of the impinging photons. Examples 
of such scattering devices are: spherical-particle suspen- 
sions (such as milk or polymer micro-spheres), polymer 
and glass multi-mode fibers and surface diffusers. Type 



The degree of entanglement and the degree of mixed- 
ness of the scattered photon pairs can be quantified by 
the tangle (T) , namely, the concurrence squared [15( , and 
the linear entropy (Sl) flq |. These quantities were calcu- 
lated from the 4x4 polarization two-photon density ma- 
trix p, by using T(/3)=(max{0, y/Ai — ^/Xz — VA3— a/Ai}) 2 , 

are the eigenvalues 
"0 -i~ 
i 



where Ai > A2 > A3 > A4 > 
of p((T2 <8> cr 2 )p* ((72 <8 cr 2 ), where cr 2 = 



and 



S L (p) = |[1 - Tr(p 2 )]. Figures [U (a) and (b) show 
experimental data reported on the linear entropy-tangle 
plane. The position of each experimental point in such 
plane has been calculated from a tomographically recon- 
structed [lH two-photon density matrix p cxp . The uni- 
form grey area corresponds to non-physical states [l~7l |. 
The dashed curve that bounds the physically admissi- 
ble region from above is generated by the so-called max- 
imally entangled mixed states (MEMS) 0, GJ]. The 
lower continuous curve is produced by the Werner states 
20] of the form: p w = VPs + ^j^Ii, (0 < p < 1), where 
I4 is the 4x4 identity matrix. Figure [2] (a) shows exper- 
imental data generated by isotropic scatterers (type I). 
Specifically, our type I scatterers consisted of the follow- 
ing categories, (i) Suspensions of milk and micro-spheres 
in distilled water, where the sample dilution was var- 
ied to obtain different points; (ii) Multi-mode glass and 
polymer fibers, where the tuning parameter exploited to 
obtain different points was the length of the fiber (cut- 
back method); (iii) Surface diffusers, where the full width 
scattering angle was used as tuning parameter. It should 
be noted that suspensions of milk and micro-spheres are 
dynamic media, where Brownian motion of the micro- 
particles induces temporal fluctuations within the detec- 
tion integration time [l4l ]. 

In Fig. 2 (a), the experimental point at the top- left 
corner (nearby T = 1, Sl = 0), is generated by the un- 
scattered initial singlet state. The net effect of scattering 
systems with increasing thickness is to shift the initial 
datum toward the bottom-right corner (T = 0, Sl = l)j 
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FIG. 2: Experimental data in the linear entropy-tangle (Sl- 
T) plane. The grey area corresponds to unphysical den- 
sity matrices. Dashed upper curve: Maximally entangled 
mixed states (MEMS), continuous lower curve: Werner states, 
(a) Polarization- isotropic scatterers (type I), (b) Birefringent 
scatterers (type II). 



the Brewster angle and (ii) Commercially available po- 
laroid sheets with manually-added surface roughness on 
its front surface to provide for wave-vector spread. All 
data thus obtained fall below the Werner curve, gener- 
ating what we called sub- Werner states, namely states 
with a lower value of tangle (T) than a Werner state, for 
a given value of the linear entropy (Sl). 




Linear Entropy (S^) 

FIG. 3: Experimental data generated by dichroic scattering 
systems (type III). 

In summary, Figs. (a)-(b) show that all data gen- 
erated by type I and II scattering systems fall on the 
Werner curve, within the experimental error; while data 
generated by scattering samples type III, which are pre- 
sented in Fig. [3J lay below the Werner curve. In Section 
III we shall present a simple theoretical interpretation for 
such results. 



that corresponds to a fully mixed state. 

Figure [2] (b) displays experimental data generated by 
birefringent scattering systems (type II). As an example 
of a system with "material birefringence" we used a pair 
of wedge depolarizers in cascade [21]. Different experi- 
mental points where obtained by varying the relative an- 
gle between the optical axis of the two wedges 22] . The 
systems with "topological birefringence" we considered 
consisted of two different devices: (i) The first one was 
a bundle of parallel optical fibers [23| • Translational in- 
variance along the fibers axes restricts the direction of the 
wave- vectors of the scattered photons in a plane orthog- 
onal to the common axis of the fibers, (ii) The second 
device was a stack of parallel microscope slides (with un- 
controlled air layers in between). This optical system is 
depolarizing because it amplifies any initial spread in the 
wave- vector of the impinging photon. This photon enters 
via a single-mode-fiber (numerical aperture NA = 0.12), 
from one side of the stack and travels in a plane parallel 
to the slides. 

In Fig. [31 experimental data generated by dichroic scat- 
tering systems (type III) are shown. We used: (i) Sur- 
face diffusers followed by a stacks of microscope slides at 



D. Error estimate 

In order to estimate the errors in our measured data, 
we numerically generated 16 Monte Carlo sets N (i = 
1, — , 16) of 10 3 simulated photon counts, correspond- 
ing to each of the 16 actual coincidence count measure- 
ments {n- xp } (i — 1, . . . , 16) required by tomographic 
analysis to reconstruct a single two-photon density ma- 
trix. Each set Ni had a Gaussian distribution centered 
around the mean value fa — n™ p , with standard devia- 
tion (ii = y/n° xp . The sets N where created by using 
the "NormalDistribution" built-in function of the pro- 
gram Mathematica 5.2. Once we generated the 16 Monte 
Carlo sets N, we reconstructed the corresponding 10 3 
density matrices using a maximum likelihood estimation 
protocol, to assure that they could represent physical 
states. Finally, from this ensemble of matrices we cal- 
culated the average tangle T av and linear entropy 5f; v . 
The error bars were estimated as the absolute distance 
between the mean quantities (av) and the measured ones 
(exp): a T = \T^ - T av |, g Sl = \S C * P - |. It should 
be noted that this procedure produces an overestimation 
of the experimental errors. In the cases where part of the 
overestimated error bars fell into the unphysical region, 
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the length of such bars was limited to the border of the 
physically allowed density matrices. 



E. Generalized Werner states 

Close inspection of the reconstructed density matrices 
generated by type II scattering systems revealed that in 
some cases the measured states represented a general- 
ized form of Werner states. These are equivalent to the 
original Werner states pw with respect to their values 
of T and Sl, but the form of their density matrices is 
different. Werner states pw of two qubits were originally 
defined [2(| as such states which are U <g> U invariant: 
pw = U <8>t7 pwU^ ®W . Here U ®U is any symmetric sep- 
arable unitary transformation acting on the two qubits. 
The generalized Werner states pgw we experimentally 
generated, can be obtained from pw by applying a lo- 
cal unitary operation V acting upon only one of the two 

qubits: paw = V <E> IpwV^ <g> I, where I = \ ? , and 



V{a,M) 



cos 7/2 



1 



sin 7/2 



("-/3) s i n7 /2 e5( Q+/3 )cos7/2 



,(1) 



In Ref. [27|, a theoretical study of the analogies 
between classical linear optics and quantum maps was 
given. Within this theoretical framework it is possible 
to build a simple phenomenological model capable of ex- 
plaining all our experimental results. To this end let us 
consider the experimental set-up represented in Fig. [T] 
The linear optical scattering element Ta inserted across 
path A can be classically represented by some Mueller 
matrix M. Q which describes its polarization-dependent 
interaction with a classical beam of light. However, Ta 
can also be represented by a linear, completely positive, 
local quantum map £ : p —> £[p], which describes the 
interaction of the scattering element with a two-photon 
light beam encoding a pair of polarization qubits. These 
qubits are, in turn, represented by a 4 x 4 density matrix 
p. Since Ta interacts with only one of the two photons, 
the map £ is said to be local and it can be written as 
£ = £a ® T, where £ a is the single-qubit (or single- 
photon) quantum map representing Ta, and X is the 
single-qubit identity map. 

It can be shown that the classical Mueller matrix M. 
and the single-qubit quantum map £a are univocally re- 
lated. Specifically, if with M. we denote the complex- 
valued Mueller matrix written in the standard basis, then 
the following decomposition holds: 



where a, f3, 7 can be identified with the three Euler an- 
gles characterizing an ordinary rotation in R 3 [24J . These 
generalized Werner states have the same values of T and 
Sl as the original pw (since a local unitary transforma- 
tion does not affect neither the degree of entanglement 
nor the degree of purity) but are no longer invariant un- 
der unitary transformations of the form U ®U . By using 
Eq. (TT]) , we calculated the average maximal fidelity of the 
measured states p C Q W with the target generalized Werner 
states p% w (p, a, We found F{p^ w ,p% w ) « 96%, 
revealing that our data are well fitted by this four- 
parameter class of generalized Werner states. 



III. THE PHENOMENOLOGICAL MODEL 



0.2 0.4 



(b) 



FIG. 4: Numerical simulation for our phenomenological 
model. Fig. 4 (a) isotropic and birefringent scattering, Fig. 
4 (b) dichroic scattering. 



T* 



(2) 



where {T^} is a set of four 2x2 Jones matrices Q, each 
representing a non-depolarizing linear optical clement in 
classical polarization optics, and {A M } are the four non- 
negative eigenvalues of the "dynamical" matrix H asso- 
ciated to M.. Given Eq. ((2|), it is possible to show that 
the two-qubit quantum map £ can be written as 



P£ = £[p] 
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fj.=0 



fj, Tfj, 



IpTt 



(3) 



where the proportionality symbol "oc" on the right hand 
side of Eq. ^ accounts for a possible renormalization to 



ensure Tr(pg) = 
essary when T4 
(i.e., dichroism). 



1. Such renormalization becomes nec- 
presents polarization-dependent losses 
We anticipate that when such re- 
normalization is necessary the map is considered non- 
trace preserving. We shall briefly discuss this issue in 
the conclusion. 

With these ingredients, a phenomenological 
polarization-scattering model can be built as fol- 
lows. First we use the polar decomposition [26j to write 
an arbitrary Mueller matrix M. = Ma-Mb-Md, where 
■Ma, -Mb and M.u represent a purely depolarizing 
element, a birefringent (or retarder) element, and a 
dichroic (or diattenuator) element, respectively. Specific 
analytical expressions for .Ma, Mb and Md can be 
found in the literature [2l[. Second, we use Eq. (O to 
find the quantum maps corresponding to A4a, -Mb and 
M d and, by using such maps, we calculate the scattered 
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two-photon state ps- In our experimental realizations 
we used isotropic scatterers Mis = Ma with isotropic 
depolarization factor < A < 1, birefringent scattering 
media M bs > described in terms of the product of 
a purely birefringent medium M b and an isotropic 
depolarizer Ma, i.e. Mbs = MbMa, and finally, 
dichroic scattering media Mds = MdMa, which are in 
turn described by a product of a purely dichroic medium 
Md and a purely depolarizing medium Ma- It should 
be noted that these product decompositions are not 
unique. Other decompositions with different orders are 
possible but the elements of each matrix might change, 
since the matrices Ma, Mb and Md do not commute. 

Filling in the above expressions with random numbers 
selected from suitably chosen ranges, we simulated all 
scattering processes occurring in our experiments. Fig. 0] 
shows a numerical simulation of the scattered states in 
the tangle vs. linear entropy plane, obtained with the 
singlet two-photon state as input state. Fig. 2] (a) 
corresponds to isotropic and birefringent scatterers, and 
Fig. 2] (b) to dichroic scatterers. The qualitative agree- 
ment between this model and the experimental results 
shown in Fig. [5] and Fig. [3] is manifest. 

IV. CONCLUSIONS 

In summary, we have presented experimental results 
on entanglement properties of scattered photon-pairs for 
three varieties of optical scattering systems. In this way 



we were able to generate two distinct types of two-photon 
mixed states; namely Werner-like and sub- Werner-like 
states. Moreover, we have introduced a simple phe- 
nomenological model based onto the analogy between 
classical polarization optics and quantum mechanics of 
qubits, that fully reproduces our experimental findings. 
In the case of sub- Werner states, the phenomenological 
model represents a non-trace preserving quantum map. 
Such description might be considered controversial since 
a non-trace preserving local map can in principle lead 
to violation of causality when it describes the evolution 
of a composite system made of two spatially separate 
subsystems (28[. However, we argue that our measured 
states do not violate the signaling condition as they 
are post-selected by the coincidence measurement, 
a procedure which involves classical communication 
between the two detectors. Finally, we expect it to be 
possible to create states above the Werner curve (in 
particular MEMS) [l8|, [l9(, by post-selective detection 
when acting on a single photon [28j. Work along this 
line is in progress in our group. 
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